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ABSTRACT 

We present the results of a search for submillimetre-luminous host galaxies of optically- 
dark gamma-ray bursts (GRBs) using the Submillimetre Common-User Bolometer Ar- 
ray (SCUBA) on the James Clerk Maxwell Telescope (JCMT). We made photometry 
measurements of the 850-/j,m flux at the location of four 'dark bursts', which are 
those with no detected optical afterglow despite rapid deep searches, and which may 
therefore be within galaxies containing substantial amounts of dust. We were unable 
to detect any individual source significantly. Our results are consistent with predic- 
tions for the host galaxy population as a whole, rather than for a subset of dusty 
hosts. This indicates that optically dark GRBs are not especially associated with 
very submillimetre-luminous galaxies and so cannot be used as reliable indicators of 
dust-enshrouded massive star-formation activity. Further observations are required to 
establish the relationship between the wider GRB host galaxy population and SCUBA 
galaxies. 

Key words: dust, extinction - gamma-rays: bursts - infrared: galaxies - stars: evo- 
lution - cosmology: observations 



1 INTRODUCTION 
1.1 Gamma-ray Bursts 

It is generally now accepted that the afterglow emis- 
sion resulting from (long-duration, soft-spectrum) gamma- 
ray bursts (GRBs) can be explained by an ultra- 
relativistic shock-wave expand i ng into a surrounding 

van Paradijs ct al. 2000[ 



medium (Meszaros & Rees 1997 



Meszaros 2001 



The precise nature of the progenitor sys- 
tems is not a settled issue. The two most popular theo- 
ries both involve stellar remnants: the collapsar/hypernova 
model, in whi ch a single mas s ive progenitor s tar undergoes 
core collapse (Woosley 1992 ; Paczyhski 199§j ) ; and the bi- 
nary merger theory, in which t wo massive stellar remnants 
such as neutron stars, merge ( Lattimer fc Schramm 1976 



Paczyhski 1986). Both of these scenarios may be able to 
explain the energetics of the explosion that produces the 
GRBs, particularly if they are beamed, since the energies of 
GRBs are comparable to those involved in the formation of 
typical stellar-progenitor black holes. 

Recently, the collapsar/hypernova model has gained 
support from three sets of observations: firstly, it was con- 
firmed that the positions of some GRBs accurately localized 
by the observation of optical and/or radio afterglows, were 
f ound to be within s tar-forming regions of their host galaxies 
(Bloom et al. 2 002), which themsel ves are frequently star- 
burst galaxies (Sokolov ct al. 2000). This is supported by 



X-ray determinations of the HI column density along the 
line of sig ht to GRBs, which is con sistent with their residing 
in GRBs (Galama fc Wijers 2001). During their creation, it 
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is likely that the massive stellar remnants required in the 
binary merger scenario would receive a substantial 'kick' ve- 
locity, so that the merger event causing the GRB would take 
place outside of the star -forming region of the host galaxy 
(Galama & Wijers 2001). Also, the delay required between 
formation of the remnants and their merger may well be long 
enough for star formation to have ceased in the host, sug- 
gesting that the host galaxy would no longer be luminous . 
Hence the position measurements of Bloom et al. (2002) 



support the single massive progenitor theory. Note though, 
that one should be wary of a potential selection effect: the 
detection of a GRB afterglow requires a certain minimum 
density in the surrounding interstellar medium (ISM), and 
so GRBs may only be identif iable when they occur inside 
galaxies (Chevalier & Li 2000). 



Secondly, in a few cases, optical afterglows of GRBs 
have been seen to contain a secondary brightening in flux 



Reichart 199S 



after a few weeks (Kulkarni et al. 199£; 



Bloom et al. 1999 



Galama et al. 200$ |Lazzati et al. 200l| ). This 



flux increase, though often not very significant, has been at- 
tributed to a supernova occurring simultaneously with the 
GRB, although alternative ideas have be en postulated, such 
as dust echoes (Esin & Blandford 200C) or interactions of 
the GRB shock wave with wind-driven density structures in 
the surrounding ISM (Ramirez- Ruiz et al. 2001). Supernovae 



associated with GRBs are compatible with the single pro- 
genitor collapsar/hypernova model rather than the merger 
hypothesis. 

Finally, the detection of iron features in the X-ray af- 
terglows suggest the presence of an iron-enriched ISM sur- 



rou nding the GRB progen itor (Amati et al. 2000 



Piro et al 



200C; Yoshida et al. 2001). This high iron mass is consistent 



with ejecta from a massive stellar progenitor. 

If the hypernova theory is the correct explanation for 
long-duration GRBs, then there should be a direct link be- 
tween GRBs and high-mass star formation activity. Since 
the gamma-ray emission from the initial explosion is not at- 
tenuated by dust, and can be detected from high redshifts, 
GRBs should be unbiased tracers pointing to star forma- 
t ion activity wherever massive stars are living and dying 
(|Krumholz et al. 199 [Blain fc Natarajan 2000|; |f3erger et al 



2001b [ frail et al. 2002| ; [Ramirez-Ruiz ct al. 2002| ) 



1.2 SCUBA galaxies 



Since the co mmissioning of the SCUBA instrument (Holland 



et al. 1999| ) in 1997 on the JCMT on Mauna Kea, a new 
era of su bmillimetre (submm) cosmology ha s been possible 
(see, e.g. Blain et al. 1999a; Smail et al. 2002] ). In particular, 
the discovery of a substantial population of dust-enshrouded 
'SCUBA galaxies' has launched a debate regarding observa- 
tional estimates of the global star formation rate in the Uni- 
verse. SCUBA traces the interstellar dust in galaxies (with 
a temperature of the order of tens of degrees Kelvin), which 
may be heated by the UV light emitted by OB stars and/or 
a hosted AGN. The relative contributions of starlight and 
AGN to dust-heating in a particular galaxy are hard to de- 
termine precisely, but studies using Chandra to observe hard 
X-rays from submm-selected galaxies suggest that ab out 20 
per cent o f the sample conta i n a detectable AGN (Bautz 
et al. 2000|; |Fabian et al. 2000|; |Almaini et al. 200l[). 



from SCUBA galaxies is powered predominantly by star for- 
mation. In this case, because of the large bolometric lumi- 
nosities of the SCUBA galaxies (at least 10 12 Lq), they make 
a substantial contribution to overall star formation activity 
in the high-redshift universe, com parable to or great er than 
that of optically-selected galaxies ( Blain et al. 1999a ) . Thus, 



they are good candidates to be hosts of many GRBs, if GRBs 
are indeed associated with the death of high-mass stars. The 
optical properties of SCUBA galaxies, where they are well- 
l ocated, indicate t hat most are very faint, with R > 25 
( [gmail et al. 20o"s ) , not dissimilar to the GRB host galaxy 
population as a whole (Djorgovski et al. 2001). It is likely 
that about 10-20 per cent of the star-formation activity 
in the high-redshift universe takes place in submm galax- 
ies brighter than 2 mjy. As a result, perhaps 1 in 5 GRBs 
should reside in such objects. This point is examined further 
in Section ^[ 

1.3 Submm observations of GRBs 

This paper describes our submm SCUBA observations of the 
host galaxies of four GRBs. To increase the chances of find- 
ing dusty SCUBA galaxies in this pilot program, we selected 
GRBs which were 'dark' in the sense that their afterglows 
were undetected in the optical despite deep, rapid searches. 
Estimates vary, but optical afterglows have been searched for 
and not found in roughly 30-50 per cent of GRBs with X-ray 
afterglows. In many cases the non-detection may simply be 
because the searches were not deep or rapid enough (Galama 
& Wijers 2001), but in some instances it is clear that the 
optical afterglows are genuinely underluminous (Groot et al. 
1998, van Paradijs et al. 2000, Lazzati et al. 2001a, Ramirez- 
Ruiz et al. 2002). An obvious possibility is that these bursts 
are heavily obscured by dust. Sin ce GRBs are expected to 
destroy any dust in their vicinity (Waxman fc Draine 200Q; 



Fruchter et al. 2001 



Reichart & Yost 2001; Venemans & 



Blain 2001), the obscuration would have to be due to dust 
elsewhere along the line of sight, as might well be expected 
in very dusty submm-bright galaxies. 

Note that we did not attempt to observe GRB after- 
glows in the submm. Target of Opportunity programs to 
observe GRB afterglow s are underway a t the JCMT and 



IRAM 30-m telescopes (Smith et al. 1999, 2001 



Frail et al 



2002). Our observations took place long enough after the 
initial explosion for only the host galaxy emission to be de- 
tectable. 

There have been two previous SCUBA detections of a 
GRB host galaxy. The first was found serend ipitously during 
a submm afterglow search (Frail et al. 2002) for a GRB that 



had an optical afterglow (GRB 010222). This suggests that 
the GRB was located in a luminous dusty galaxy, but either 
in a relatively dust-free region or near the edge of the galaxy, 
so that the ex plosion could clear ou t any dust along the line 
of sight. Also, Berger et al. (2001b) reported the results of 



a successful targeted search for submm emission from GRB 
000418. This too had an optical transient identification. An- 
other relevant detection was that of a GRB host in the radio 



(GRB 980703; Berger et al. 2001a), whose flux was claimed 



Therefore, it is likely that the luminous dust emission 



to be caused by a large star-formation rate rather than AGN 
activity, based on optical spectroscopy and the absence of 
radio variability. However, a deeply embedded AGN within 
a dust- and gas-rich galaxy could plausibly contribute the 
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Table 1. Photometry observations log. See notes about GRB 
001025 and GRB 001109 in section §. 



Source name 


Date of obs. 


Integ. time (s) 


T 850/im 


GRB 970828 


24/04/01 


900 


0.23 


GRB 981226 


03/10/01 


2052 


0.27 




24/11/01 


2115 


0.27 


GRB 990506 


08/05/01 


2250 


0.26 


GRB 000210 


01/09/01 


1800 


0.24 




03/10/01 


2250 


0.27 


GRB 001025 


05/10/01 


2250 


0.30 


GRB 001109 


08/05/01 


1350 


0.26 




03/10/01 


2700 


0.26 



emission from this galaxy, especially as the radio source is 
located very near to the nucleus of the host galaxy. 



2 OBSERVATIONS 

Observations were mad e using the 850-^tm photometry pixel 
on the SCUBA array ( [Holland et al. 1999| ). They are sum- 
marized in Table [j]. Simultaneous observations were made 
with the 450-/im photometry pixel, but weather conditions 
and/or array noise were never good enough to yield use- 
ful data. Observations were made in the standard photom- 
etry mode, using a 7 Hz, 60" 'chop' in azimuth to provide 
a blank-sky reference and a further telescope 'nod' to pro- 
duce a measure of any sky gradient. Each set of photometry 
observations takes 18 seconds. Column 3 of Table 1 shows 
the integration time on-source for each GRB. Pointing was 
checked regularly and was always better than 2.5", much 
smaller than the primary beam of the array of 14" at 850 /im. 
Seeing was monitored, especially around dawn and dusk, and 
data was only taken when the seeing was less than 1". Sky 
transparency was calculated by interpolating between regu- 
lar skydip values, though the situation was monitored more 
frequently with data from the JCMT's water vapour moni- 
tor. The average sky opacity at 850 fim during each source's 
observations, r sso^m, is given in column 4 of table 1. For flux 
calibration, identical observations were made of the planets 
Mars and Uranus, and where necessary the secondary flux 
calibrator CRL 618. We found no significant deviations be- 
tween the observed and expected fluxes of these standard 
objects. 

Unfortunately, since the observations were made, two 
of the GRBs in our sample have been found to have been 
wrongly located: GRBs 001109 and 001025. In the case of 
GRB 001109, the radio-located position which we observed 
i s now thought to corr espond to a faint constant radio source 
(Berger & Frail 2001), and thus to be a mis- identification 
of the GRB afterglow. We detected a net positive flux of 



(1.89 ± 1.40) mjy from this source, consistent with it s likely 
identification as a high-redshift star-forming galaxy ( |Berger 
fc Frail 200 For GRB 001025, XMM error box SI was 



observed, since it contained a candidate host galaxy (Hjorth 
priv. comm.), but subsequently Hurley (priv. comm.) has 
calculated that this error box lies outside a revised IPN an- 
nulus. Hence the results of the observations of both these 
objects are not included in the rest of this paper. 



Table 2. Coadded 850-fim results for the four objects observed. 
Redshifts z, and R magnitudes, where determined, for candi- 
date host galaxies are shown here ( Djorgovski et al. 2001 (GRB 
970828), |Holland et al. 200o| CGRB 981226). bloom ct~al. 2001 



(GRB 990506) and [Piro et al. 2002| (GRB 000210)) 



Source 


Flux (mjy) 


z 


R mag 


GRB 970828 


1.26 ± 2.36 


0.96 


25.5 


GRB 981226 


-2.79 ± 1.17 




24.3 


GRB 990506 


-0.25 ± 1.36 


1.3 


24.8 


GRB 000210 


3.31 ± 1.54 


0.85 


23.5 



3 DATA REDUCTION 

The data w ere reduced using the stan dard Starlink SURF 
procedures (Jenness & Lightfoot 200C). Particular care was 
taken over the removal of atmospheric noise, since the ex- 
pected low fluxes of the sources make this an important 
factor. This sky-noise removal is possible using the other 
bolometers on the array, as suming that they are pointing 

the median 



tsaak et al. 



(2002p , 



at the blank sky. Following 
value of the signal from the reference bolometers was used 
rather than the mean. Usually the reference bolometers used 
are the inner ring on the 850-/im array. For the October 
2001 observations (see table ^), however, the whole array 
suffered from elevated noise. Particular care was taken to 
choose only those bolometers with normal noise levels for 
this final sky-removal stage, and the bolometer noise values 
were measured far more frequently than usual. The results 
presented have been clipped at the 3<r level, though this 
generally had little effect on the final results. 



4 RESULTS AND DISCUSSION 

Table ^| shows the overall weighted mean at 850 fj,m for each 
of the four reliable sources observed. The overall, weighted 
mean flux for all four sources is -0.37 ± 0.82 mjy. 

Models for the evolution of the star formation rate in 
dusty galaxies can be used to pre dic t the likely 850- Atm fluxes 
of GRB hosts (|Blain et al 1999cj ,|b|; |Blain fc Natarajan 200c| ; 
Ramirez- Ruiz et al. 2002), as shown in Fig. [j]. We assume 
that: dust in SCUBA galaxies is predominantly heated by 
high-mass stars; GRB rates are tied to the rate of formation 
of high- mass stars; and most high-redshift star formation ac- 
tivity is enshrouded by dust. Results for two models based 
on infrared and submm data are presented, a simple para- 
metric model of th e evolution of low-redshift galaxies (BSIK: 
Blain et al. 1999a ) and a model ba sed on luminous h icrarchi- 



cal merging of galaxies (BJSLKI: Blain et al. 1999b). These 
models both provide a good representation of the evolution 
of dusty galaxies, measured in other observations. About 20 
per cent of all GRB host galaxies are expected to have fluxes 
above the SCUBA confusion limit for detection of 2mJy; 
about 10 per cent are expected to have fluxes greater than 
5mJy. In general, we predict an average across all sources 
of 0.9 mjy. The solid histogram shows the result of our ob- 
servations, and the dotted histograms define the ±1<t errors 
in flux. 

It is clear from Fig. |l| that the errors on our results 
are consistent with the predictions of the galaxy evolution 
models, and that the overall mean flux is within 2a of the 
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1 2 
Host flux density S / mJy 



Figure 1. Compa rison between results and predictions ( Ramirez- 
Ruiz et al. 2002). The smooth curves show the fraction of all 
GRB host galaxies expected to exceed a given 850-^tm flux i n two 
current mod els — for a pure luminosity evolution (BSIK: 31ain 



et al. II )99a ) and an hierarchical model respectively (BJSLKI: 
Blain ctah 1999kJ). Stepped lines show results from this pilot 



tical and radio afterglow behaviour may all be explained by 
a complicated density structure in the ISM around the GRB. 
In particular, a cavity in the ISM density would explain the 
rapid X-ray and radio decays, and would predict a simi- 
larly fast optical decay, meaning that the optical afterglow 
searches were too slow/shallow to detect the emission. Such 
a structure could be produced by the mass loss phases that 
massive stars (which are of course plausible GRB progeni- 
tors) are known to go through in their post main s equence 
lives ( Chevalier fc Li 200C ; Ramirez-Ruiz et al. 2001 ). Hence, 
there may be no need to infer the presence of dust around 
this GRB. Other suggestions for rapid decays have been pro- 
posed which also do not rely on the presence o f dust, su ch 



as the effects o f jet structure (Frail et al. 2001 



fc Kumar 2001) 



Panaitescu 



Our non-detection is consistent with the brief weather- 



affected SC UBA afterglow search for this GRB by Smith 



study - see text. 



4.1.2 GRB 990506 

GRB 990506 appears to have some similarities with GRB 
981226. Again it ha d a very rapidly decaying radio transient 
( Taylor et al. 200C ) . T he optical transient was not detected 
at R = 19 after 1 hour (|Zhu fc Zhang 199S| ), nor to R = 23.5 
after 11 hours (Pedersen et al. 199£ ) . The rapid decay of the 
radio afterglow again suggests that a non-dusty effect, local 
to the burst, may explain the absence of an optical afterglow. 
Optical searches for t he host have identif ied it as a very faint 
and compact galaxy (Bloom et al. 2002). 



expected value. GRB 000210 is clearly our best candidate for 
a detection, but only at a significance of about 2a. However, 
observations of GRB 000210 suffered from the raised array 
noise discussed in section [], and from its low elevation from 
Mauna Kea. Hence, the atmospheric noise for GRB 000210 
was greater than for our other sources. 

Note that we purposefully chose a selection bias that 
should have increased our chances of finding dusty GRB 
hosts, by observing 'dark bursts'. Since we did not see any 
increase in our detection rate, alternative explanations are 
necessary for the lack of optical afterglow detections. We 
now look briefly at the four reliable GRBs in our sample in 
more depth. 



4.1 Individual Sources 

4.1.1 GRB 981226 
Three e arly candidates for the optical afterglow of GRB 



4.1.3 GRB 970828 



981226 (|Castro-Tirado et al. 199$ |Galama et al. 1998; ; |Woz- 
niak et al. 1998|) were reject ed upon the location of a radio 
afterglow ( |Frail et al. 1999 ). The limit on the optical af- 
t erglow is therefore R > 23.5 at 10 hours after the alert 
([Lindgren et al. 1998j). Later HST/ST IS and VLA imagin g 



has located the probable host galaxy (Holland et al. 200C) 



The multi-wavelength afterglow emission from GRB 
981226 provides some clues. The X-ray afterglow was found 
t o have a double-pea ked str ucture followed by a rapid decay 
( |Frontera et al. 2000[ ). Also, |Frail et al. (1999 ) noted a rapid 
decline in the radio afterglow. Taken together, the X-ray, op- 



GRB 9708 2 8 was dark to a depth of R = 23.8 ( Groot 



ct al. 199£ ; |Djorgovski et al. 2001 ) in observations taken 



from 4hrs after the burst. A ra dio flash was observed by 
the VLA (Djorgovski et al. 2001) and has been interpreted 
as reverse shock emission. No conventional radio afterglow 
was detected. Subsequent Keck and HST observations of 
the location of th e radio flare revealed an interacting three- 
component host (Djorgovski et al. 2001), with a possible 



identification of the GRB location in a dust lane between 
the two brightest components. In calculations considering 
both jet and spherical models for the GRB shock geome- 
try, Djorgovski et al. (2001) conclude, based on the X-ray 
afterglow flux, that a single typical Giant Molecular Cloud 
could provide all the extinction necessary to fit the upper 
limits to the optical afterglow flux density. They also note 
that the Keck and HST images indicate that the two host 
galaxy components on either side of the GRB location are 
both slightly but not highly reddened, suggesting a low total 
dust mass in the system.rl Our positive but not significant 
measurement supports the hypothesis that in this case at 



There are however indications from observations of galaxies 



know n to be very luminous at submm wavelengths (Ivison et al 



2001) that regions where dust emission is strong may not correlate 



with regions which at optical wavelengths appear to be reddened 
by strong extinction. Hence the lack of obvious extinction in op- 
tical images does not rule out the presence of a large amount of 
illuminated and heated dust. 
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least, it was not galaxy- wide dust that caused the obscu- 
ration of the optical afterglow, but rather a localized cloud 
or clouds of dust along the line-of-sight, consistent with the 
X-ray results. Of course, SCUBA's resolution is such that 
our observations are sensitive only to the total (illuminated) 
dust in the entire system. 



4.1.4 GRB 000210 

GRB 000210 represents our most likely host galaxy detec- 
tion. Both X-ray and radio transients were found for this 
source, but no optical transient (R > 23.5 at 12.4 ho urs af- 
ter the burst, (|Gorosabel et al. 2002|; |Piro et al. 2002|)). The 



X-ray transient did not display the rapid decays found for 
GRBs 981226 and 990506, leading Piro et al. (2002) to reject 
the no-dust hypotheses discussed above. Instead they con- 
clude that the most likely scenarios are either obscuration 
of the optical transient by a clumpy local environment, or 
line-of-sight obscuration by the whole host galaxy, either of 
which is allowed by our findings. 



5 CONCLUSIONS 

From our small sample of four reliably-identified dark GRBs, 
we find that the 'dark bursts' do not preferentially select 
dusty host galaxies with very significant amounts of star 
formation. Looking at some members of our sample, we can 
explain the lack of optical transients for other reasons. It 
seems likely that to characterise the optically dark bursts as 
a physically distinct population of GRBs would be mislead- 
ing. In each case different circumstances due to a combina- 
tion of the observing conditions and the physical conditions 
at the location of the GRB could give rise to the lack of an 
optical afterglow. If, instead, we view the GRBs included 
here as four examples of the overall population, we find that 
the results agree with the distribution predicted assuming 
that GRBs tracew high-mass star formation. Our lack of a 
single strong detection implies that no more than 20 per 
cent of GRB hosts are submm galaxies detectable at a flux 
density brighter than 5mJy using SCUBA. In interpreting 
the results, however, the small sample size inevitably makes 
it hard to draw any solid conclusions. 

Combining our three most sig n ifican t results with the 
observations of Smith et al. (1999, |200lD , we find that 11 
GRB locations have been observed with SCUBA to noise 



levels < 2m Jy, and in all but two (GRB 980329, Smith et al 



199E and GRB 000210, this paper) detection of submm emis- 



sion from a host gala xy was ruled out . The two GRB host s 



detected by SCUBA (Frail et al. 2002 Berger et al. 2001b) 



had afterglows with optical transients. Therefore, while some 
GRBs definitely are located in dusty galaxies, the route to 
selecting these GRBs on the basis of their afterglow data is 
not yet clear. Sin ce GRBs can only clea r dust out to less than 
100 pc distance (Fruchter et al. 2001), it would be surpris- 



ing if optically-selected samples of GRBs are not generally 
biased against dusty hosts. 

However, as noted above, there may be alternative ex- 
planations for the optical faintness of some afterglows be- 
yond location in a dusty host galaxy. Our results indicate 
that radio-located optically-dark bursts seem not to be re- 
liable indicators of luminous, dusty host galaxies. It may 



be that the physical conditions of the ISM in the densest 
star forming regions are incompatible with the generation 
of intense radio emission from GRB shocks. In that regard, 
it would ultimately be interesting, and certainly possible in 
the SWIFT era, to study a sample of hosts of purely X- 
ray selected GRBs, with the hope that the prompter X-ray 
emission would be less affected by the wider environment of 
the progenitor. 

From the SCUBA observations alone, we cannot sepa- 
rate the possibility that there is little dust in these systems 
from the possibility that there is a lot of dust but insuffi- 
cient UV photons to make it glow brightly. Alternatively, 
the dust may be heated to high temperatures and so can- 
not be detected by SCUBA at all. The predicted flux of 
a 5 x 10 12 Lq galaxy at z = 1 decreases from 20mJy to 
0.08 mjy as the dust temperature increases from 20 K to 
80 K (see Fig. 5 in |Blain et al. (2002] )). Typical dust temper- 
aturcs for submm -selected objects are thought to be 40 K 
(fvison et al. 2000), but a temperature ~ 60 K would be hot 
enough to prevent such a galaxy being detected above the 
2mJy confusion limit of SCUBA. However, optical colours 
of the host galaxies of GRBs 970828, 9 8 1226 and 000210 
(|Djorgovski et al. 200l|; |Frail et al. 1999|; |Piro et al. 2002j) 



show only modest reddening, to the optical depths that can 
be probed, and thus perhaps imply little dust overall in the 
hosts. Ideally, high-resolution observations to locate dust- 
enshrouded star-formation activity with respect to the loca- 
tion of GRBs need to be made. 

The results presented here suggest that a more extensive 
deep submm survey of GRB hosts is necessary to investigate 
their far-infrared and submm properties in detail, and this 
is now underway at the JCMT (Barnard et al., in prep.). 
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